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<— I ■ This talk presents details of the calculation of two-loop corrections to heavy 

Qh| baryon currents containing one heavy quark in leading order of the l/mn-expansion, 

i.e. in the limit of Heavy Quark Symmetry (HQS). The calculations lead to the 

determination of the anomalous dimension of these currents. I also discuss problems 

with different 75-schemes and their solution by the finite renormalization procedure. 
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1 Introduction 

Heavy baryons such as A c , E c , A^ or E& containing one heavy quark can be successfully 
described by the Heavy Quark Effective Theory (HQET), especially in leading order of 
the 1/mQ-expansion, i.e. in the limit of Heavy Quark Symmetry (HQS) (for a review on 
HQET see [|l|, for details on HQS see [[|). For each of the ground-state baryons there are 
two independent current components, e.g.Q 

J ax = [q iT CT-f 5 q j }Q k e ljk and J A2 = [q iT CT^ q j }Q k e ijk , (1) 

while the general structure of heavy baryon currents has the form 

J = [ q ^CrT^]T'Q k e ijk , (2) 

where the index T means transposition, C is the charge conjugation matrix, r is a matrix 
in flavour space, and i, j and k are colour indices. The effective static field of the heavy 
quark is denoted by Q, and T and V are the Dirac structures of the different current 
parts. 

In Fig. 1 are shown some examples of diagrams up to two-loop order which we calcu- 
lated. The one- and two-loop diagrams can be grouped in three different classes, namely 
the heavy diquark system (diagrams (b), (c), (e), (f) and (g)), the light diquark system 
(diagrams (d) and (h)) and the irreducible vertex contribution (diagrams (i) and (j)). The 
first two can be handled in a manner similar to the meson case, and actually it turns out 
that the replacement Cb —>■ C F leads from the diquark to the meson quark-antiquark case 
(and not vice- versa). So the mesonic case can be seen as a byproduct of our calculations. 

2 Algorithmic calculations 

The calculation of all these diagrams can be executed by nearly the same procedure. 
The procedure is strictly algorithmic, so that the calculations could be automatized and 
programmed in MATHEMATICA. In the following I want to emphasize the main steps of 
this calculation, and I want to do this for the first two-loop heavy diquark diagram which 
is displayed in Fig. 2 (the arrows on the gluon lines only indicate the momentum flow). 

2.1 Colour structure 

Before starting the main calculation, the colour structure can be treated. In general it 
contains the Gell-Mann matrices, the epsilon tensor for resons of colourless, and the colour 
couplings of the gluons, which consists of a 5 a b for a single gluon and a factor ij a bc for a 
three-gluon vertex. In our case we get 

qV\T%(T%Q k '' (T a %(T%e ijk = (^^) 2 qWQ k e ijk =: Clq^Q k e ljk . (3) 

N c is the number of colours, and the results can be expressed by the constants Ca = N c , 
C B = (N c + l)/2N c , C F = {N% - l)/2N c , T F = 1/2 and N F , which denotes the number 
of flavours. 
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2.2 Other structures 

Without the colour structure we get 
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The four-momentum of the heavy quark is given by q = itiqv+p, and we have u = pv . This 
vertex correction can be split up into a Dirac structure Tl)^f ll ili r ) u T 1 a momentum structure 
kH v and a scalar integral (also called "race integral") 

tihi) gt f ff-l\ 2 f-l\( -1 \ w uj al n k d n l 
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2.3 One- and two-loop functions 

With the help of the recursion formulas developed in 0, [|, || it is possible to reduce 
each two-loop function 7 5 of the HQET and each two-loop function G 5 of the massless 
QCD to corresponding one-loop functions I<i and G2, which dissolve into Euler's Gamma 
functions. The recursion formulas are one of the main parts of the developed program 
package, and we would like to thank David Broadhurst for providing us with a REDUCE 
implementation of these formulas. 

2.4 Momentum structure, covariant evaluation and operators 

By extracting the Dirac structure, the remaining integral receives corresponding Lorentz 
indices. It can be evaluated covariantly, in this case into 

AgV + Bv»v v . (6) 

The parts A and B can be obtained by contracting the integral with the corresponding 
dual tensors. This contraction pairs the momentum structure to the scalar products 

M = -(k 2 + l 2 -(k-l) 2 ) and (kv)(lv), resp. (7) 

The scalar products "operate" on the entries of the two-loop function (e.g., k 2 decreases 
the first entry by one unit), and can therefore be identified with operators on the different 
entries of the two-loop functions. 



2.5 Dirac structure and its reduction 

The covariant evaluation in combination with the Dirac structure gives 

b[ hl) = i, lph i, lv Y{Ag^ + Bv^v") = At-i^T + BfttfT = ((2 - n)A + B)Y. (8) 

For the heavy diquark system the matrix V will always remain on the edge of the structure. 
In contrast to this, for the light diquark system the Dirac structure can be reduced to one 
of the main structures 

r := r, r\ := 7 ^rw\ r 2 := 7 ^r 7 v, (9) 

r 3 := l i avl P 't^i'l P l v l il and T 4 := ~i v n v ~i p ~i (T Tl' J l P l v l 1 ' '■ 

3 Result of the perturbation series 

The result of the perturbation series to second order for the different ground-state hadron 
currents can be written in quite a compact form, if we make some slight assumption for 
T. In this case, we can split off a scalar valued vertex junction V. 

3.1 Introducing n 7 and s 7 

If T is totally antisymmetrized, T e {1,7^,7^7^, . . .}, we get 

7 M IV = (-ir(n-2n 7 )r and 7o r 7 o = (-l) n ^ 7 r, (10) 

where n is the space-time dimension, n 7 the number of Dirac matrices and s 7 denotes the 
number of 70's in T, i.e. being +1 for an even and —1 for an odd number of 70's. Problems 
will emerge from the appearence of 75, but I will postpone the discussion of this to the 
end of the talk. 

3.2 Vertex function for the heavy baryon currents 

Without looking at the parts from different classes of diagrams, the vertex function for 
the heavy baryon current in its generalized form is given by 

V = 1 + 2AV {hl) + AV {ll) + AV {ir) (11) 

This expression can be evaluated in powers of 1/e with e = (4 — n)/2 and shows singular- 
ities which can be canceled by means of the minimal subtraction scheme (MS). Indeed, 
the evaluation in 

- i :=-- lE + ln(47r) = - ln(47re 7B ) (12) 

is rather more simple. This evaluation allows for the application of the modified minimal 
subtraction scheme (MS). By using a general gauge parameter a (a = 1 means Feynman 
gauge) for the one-loop contributions, we get 

00 fa \ m °° m (a \ m 1 

"=E(£) v-=EE(S) h v ~ wlth (13) 

m=0 \<m/ m=0k=0 \^n/ t 



V; 1 = C B ((n 7 -2) 2 + 3a-l), 

V? = C B (2n* -6n 7 + 3 + a+(n 7 -2)s 7; . 

V 2 2 = ^(3C B ((n 7 - 2) 2 + 2) 2 - C A (lln 7 - 44n 7 + 51) + 

+4JV F T F (ra 7 - 3)(n 7 - 1)) und (14) 



V2 = ^(9C B ((n 7 - 2)(13r£ - 70n 2 + 126n 7 - 60 + 4(n 2 - 4n 7 + 7)s 7 ) + 32((2)) + 
36 v 

-C A (18n 7 - 144n 7 + 289n 2 + 128n 7 - 516 + 72((2)) + 
-72C F (n 7 - 3)(n 7 - 1) - 4iV F T F (n 2 - 16n 7 + 24)). 

V 2 ° is omitted, because it has no significance for the following calculations. 

4 The renormalization 

The vertex function still contains UV- singularities. It is a bare quantity, from now on 
denoted by V°, which has to be renormalized. This renormalization is done with the help 
of a renormalization factor. Dimensional regularization is used from the beginning, what 
remains is the subtraction with the aid of the modified minimal subtraction scheme, as 
indicated above. In the following we simplify the notation by using the term e instead of 
e', but we keep in mind that we have used the MS-scheme. 

4.1 From the perturbation series to the renormalization factor 

The divergencies of V° will be absorbed by the renormalization factor Zy. We have to 
state therefore that Z^V = V is a finite quantity. By expanding the renormalization 
factor Zy in a double series 

Zv = l+EE(^) :i^ = E-> fc with Z° = l, (15) 

m=lk=l v ^" 7 fc fc=0 fc 

a comparison of the coefficients up to second order gives 

Z\ = V ± , Z\ = K> 2 and Z\ = Vj - V?V?. (16) 

4.2 Renormalization of coupling and gauge parameter 

Up to now we have left out the renormalization of the coupling constant oc s and the gauge 
parameter a, 

a° s = a. (l + ^S a + 0(a 2 )) , a = a (l + ^-5 a + 0(a 2 )) . (17) 

Looking at terms up to second order in perturbation theory, we only need to correct 
the terms of first order. In this way we get corrections for the coefficient Z\ of the 
renormalization factor Zy, namely 

Z\ = V, 1 - VfV, 1 + S a V° + 8 a V?, (18) 

where V x is the gauge dependent part of V® . 



4.3 From the renormalization factor to the anomalous dimen- 
sion 

The renormalized coupling is an object in n space-time dimensions. Therefore, it gets a 
double dependence in the renormalization scale /x, 

a s = a s (fx,e) = a s (/i)/i 2£ (a s (/i) is the coupling in 4 space-time dimensions). (19) 

The beta-function for the coupling and the gauge, two essential tools for the renormal- 
ization group equation, can then be calculated as 

_ dlna 8 (n,£) _ d\na s (fi) dlnajfj.) 

a In /j, dln/j dux ft 

This can be used for the calculation of the anomalous dimension 7, 

dZ((Xg(fi),a(iJL);e) dZ dhia s (/i) dZ dlna(/i) 



1 Z 



dhi/i dlna s dha.fi d In a din /x 

" Z (/3-2.) + #-A. (21) 
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By evaluating the finite quantities (3, (3 a and 7 in a perturbation series, starting at first 
order, the comparison of the coefficients leads to a determining equation for the coefficients 
of the anomalous dimension, 

7m = -2mZ l m , (22) 

and on the other hand to a class of consistency checks, which can be revolved to the 
coefficients of the perturbation series, e.g. 

2V 2 2 = V'Vl + 8 a Vl + 5 a V 1 1 (23) 

where VJ 1 again denotes the gauge dependent part of V±. For the derivation we should 
note that 

8 a = \Px = ~C A + \N F T F and 5 a = U(5 a ) x = ^--^C A - ^N F T F . (24) 
2 6 6 2 06 

5 The anomalous dimension of the 
generalized heavy baryon current 

To renormalize the current, we also have to include the renormalization of the outer legs, 
because we have 

jo = [ q 0T CrV Tq }T'Q = Z q Z% 2 Z v J = ZjJ => lJ = 2 lq + lQ + lv . (25) 

By taking this into account || |7|], to first order we get 

7j = | i ( - 2C B ((n 7 - 2) 2 + 3a - 1) + 6C F (a - 1)) + 0(« 2 ). (26) 

For the values Cb = 2/3 and C F = 4/3 of SU(3)c, this expression is gauge independent. 
So we can use these values, C A = 6 and T F = 1/2 and restrict ourselves to the Feynman 
gauge. 
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5.1 Different 75-schemes 

The problem in presenting the final result consists in the occurence of different schemes 
for 75. 



In the scheme due to 
lohner and Maison [Ft 



't Hooft, Veltman, Breiten- 
B], an additional 75 in the 



vertex chances n 7 to 4 — n 7 and s 7 to — s 7 . For the 
naively anticommuting 75-scheme [1(| , both defin- 
ing equations for n 7 and s 7 change their general 
sign, but this will cancel for all main structures IV 
The values for the parameters for different heavy 
baryons are displayed in the enclosed table. 
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s 7 


baryon 


7 5 A ° 
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7 5 AC 7o 


1 


-1 


A 2 


7 


1 


+1 


Si,S* 
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2 


-1 


S 2 , S2 


7 5 HV 


4 


-1 


Ai 


7 5 HV 7o 


3 


+1 


A 2 



5.2 Values for the anomalous dimension 

For the anomalous dimension of different baryon currents and different schemes we get 



7ai 



7$ 



7si 



7S2 



-/ HV 

7A1 



7% 



^)+i(16C(2)+40JV F -796)(^) : 

- 4 (£) + 5 (16C(2)+20iV ''- 322) (S) ; 
- 4 (£) + 5 (16C(2)+20Ar ''- 290) (S) ; 
-5(£) + >^) +8iV ^ 324 )(£ 

^)+i(16C(2)-24 J V F + 260)(^) ; 

- 4 (£) + 5 (16C(2) - 12JVf+206) (S) ; 



(27) 



The "naive non-Abelianization" recipe, the estimate of the total result by calculating the 
quark-loop diagram and replacing Np by Np — 33/2, does work for most of the cases with 
an accuracy of 15 to 30%. Only in the case of the current J^ 2 it seems to fail. Note, 
however, that in this case the anomalous dimension is one order of magnitude smaller 
than in all other cases. 



5.3 The finite renormalization 

The difference of the anomalous dimension for the two different 75-schemes turns out to 
be a correction of second order which is proportional to the coefficient (3\ of the beta- 
function of the coupling. So we might think about connecting the already renormalized 
currents by some finite renormalization factor |T1|, namely 



HV/ 



j ao Gu) = z r j H ». 
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So we get 
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(29) 



Indeed, the ansatz Z? = Ar(a s (/i)/47r) leads to Arfii = Or with 

16 „ 



3 ~ 

(30) 



"^75 


— — w^xxw^ t«i< , -ij<7, 


^75 7i = ^7570 


= -C B (llCA-4iV F r F ) and 


^Iblilj = ^757i70 


= o, 


and therefore 




7—1 8 f° 7 — 7 

^75 — - 1 - ° a u Bi ^757i ^7570 


Ot 

= 1 — 4—C B and Z^ 7ilj = Z 757i70 



0. (31) 

The vertex structures T = 7571, 757i7j and 7s7i7o are included for the sake of complete- 
ness, although they are not needed for the baryonic s-wave states discussed in this talk. 
The meson case, already treated by Trueman and Larin [II |, can again be obtained by 
the replacement Cb — > Cp. 
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Figure Captions 



Fig. 1: Examples for corrections up to two-loop order to the baryonic vertex. Dou- 
ble lines indicate heavy quarks, single lines indicate massless quarks, the 
springs stand for gluons 

Fig. 2: The first heavy diquark diagram (Fig. 1(e)) with momentum labels 
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